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I  .  INTRODUCTION 


A .  PURPOSE 

This  report  describes  an  experiment  intended  to  detect  possible 
electron-density  perturbations  in  the  ionosphere  associated  with  small 
travelling  fluctuations  of  the  earth's  magnetic  field.  The  estimated 
size  of  the  perturbations  and  the  problems  involved  in  their  detection 
are  discussed. 

B .  BACKGROUND 

In  1722  George  Graham,  a  London  instrument  maker,  observed  that  a 
compass  needle  is  never  completely  at  rest  but  instead  oscillates 
continually  about  its  mean  heading.  It  was  thus  discovered  that  the 
earth's  magnetic  field  is  in  a  continual  state  of  flux.  The  nature  of 
this  variability  has  been  the  subject  of  considerable  investigation  for 
the  ensuing  two  and  one-half  centuries. 

Many  observed  magnetic  fluctuations  have  a  quasi-periodic  structure. 
These  form  a  group  of  disturbances  called  geomagnetic  micropulsations, 
the  "micro"  prefix  stemming  from  the  fact  that  their  amplitudes  are  very 
small  in  relation  to  the  mean  value  of  the  earth's  field.  A  further 
classification  sorts  these  micropulsations  into  a  number  of  subgroups 
based  on  their  morphological  properties.  This  report  is  concerned  with 
particular  properties  of  the  Pc  1  micropulsations,  a  subgroup  composed 
of  regular  magnetic-field  oscillations  having  frequencies  in  the  range 
0.2  to  5  Hz.  At  middle  latitudes  these  oscillations  are  primarily 
nighttime  phenomena  having  amplitudes  ranging  from  a  few  milligammas  to 
hundreds  of  milligammas,  and  persisting  from  a  few  minutes  to  several 
hours.  The  classification  of  micropulsations  and  the  units  in  which 
they  are  measured  are  discussed  in  more  detail  in  Appendix  A. 

Pc  1  micropulsations  apparently  originate  as  plasma  instabilities 
in  the  equatorial  magnetosphere,  from  which  region  they  travel  down  to 
the  earth's  surface  along  high-latitude  magnetic  field  lines  as 
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hydromagnetic  waves..  Figure  i  is  .  schematic  representation  of  this 
situation.  Given  no  other  propagation  mechanisms,  it  might  be  expected 
that  the  micropulsations  would  only  be  observed  near  the  points  where 
the  disturbed  magnetic  field  lines  penetrate  the  earth's  surface.  How- 

ever,  Pc  l's  are  observed  worldwide  e-u-  * 

° iawiae ,  although  the  strongest  and  most 

frequent  occurrences  are  indeed  at  high  latitudes. 

Of  particular  interest  is  the  observation  that  the  occurrence  of  a 
Pc  1  micropulsation  outburst  at  a  middle-latitude  location  is  often 
correlated  with  similar  activity  at  high  latitudes  on  the  same  geo¬ 
magnetic  meridian,  in  addition,  the  middle-latitude  event  is  usually 
observed  to  be  delayed  in  time  with  respect  to  the  corresponding  high- 
latitude  event,  as  if  the  energy  causing  the  disturbances  had  travelled 
With  a  measurable  velocity  from  north  to  south  (in  the  northern  hemis¬ 
phere)  .  The  propagation  velocity  which  results  from  a  consideration 
of  the  time  and  distance  between  the  observations  turns  out  to  be 
approximately  the  group  velocity  of  a  hydromagnetic  wave  of  Pc  1 
frequency  travelling  i„  the  ionosphere  near  the  region  of  maximum 
electron  density.  Hence  the  theory  has  arisen  that  Pc  1  micropuls. tlon 
energy  travels  fro.  high  to  middle  latitudes  as  a  hydromagnetic  wave 
and  that  a  sig„ificant  portion  of  the  wave  energy  is  confined  to  the 
ionosphere.  This  concept  is  also  shown  schematically  in  Figure  1 

A  more  detailed  discussion  of  the  theory  describing  the  creation 
and  propagation  of  micropuls.tions ,  as  well  as  further  information  on 
their  properties  and  classification,  are  given  in  Jacobs  (1970),  He 
reviews  the  whole  micropulsation  field  and  includes  a  good  bibliography 
in  addition,  a  good  account  of  the  early  history  of  geomagnetism  can  be 
found  in  Chapman  and  Bartels  (1940) 


METHOD  OF  OPERATION 


A  hydromagnetic  wave  with  a  frequency  of  about  one  Hertz  travellin. 
through  a  plasm,  with  the  attributes  of  the  earth's  ionosphere  should 
indue,  in  this  Plasma  a  periodic  electron-density  perturbation  having 


* 


Some  of  the  to --.is 
wave,  ionosphere 


' ntroduced  in  this  section  (plasma,  hydromagnetic 
and  magnetosphere)  are  defined  in  Appendix  B. 
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the  same  frequency  (Fraser-Smith ,  1969).  If  such  a  perturbation  could 
be  observed  in  the  ionosphere  while  micropulsations  were  being  observed 
on  the  ground,  then  support  would  be  lent  to  the  theory  that  the  micro¬ 
pulsation  energy  flows  from  high  to  middle  latitudes  primarily  in  the 
ionosphere. 

A  high-frequency  (HF)  radio  wave  launched  vertically  from  the 
ground  will  be  reflected  by  the  free  electrons  of  the  ionospheric  plasma 
(provided  that  the  wave  frequency  is  not  higher  than  the  critical 
frequency)  and  will  return  to  the  vicinity  of  the  transmitter.  At  the 
ground,  the  phase  of  the  reflected  wave  will  vary  with  respect  to  the 
phase  of  the  transmitted  wave  if  the  electron  density  along  the  trans¬ 
mission  path  varies  with  time.  In  particular,  a  periodic  electron- 
density  perturbation  induced  by  the  passage  of  micropulsation  energy 
should  produce  a  periodic  phase  modulation  of  the  reflected  wave. 
Detection  of  periodic  electron-density  perturbations  in  the  ionosphere 
is  thus  reduced  to  detection  of  phase  modulation  of  a  radio  wave. 

An  experimental  system  designed  to  detect  Pc  1  micropulsation 
energy  in  the  ionosphere  in  the  manner  described  ebove  was  set  up  by  the 
Stanford  University  Radioscience  Laboratory  in  mid-1970.  A  low-power 
HF  wave  was  transmitted  from  a  site  near  Pescadero,  California,  and  was 
received  on  the  university  campus  a  few  miles  away  after  being  reflected 
by  the  ionosphere.  The  theoretical  and  practical  aspects  of  this  project, 
as  well  as  the  experimental  results,  are  discussed  in  the  following 
pages. 
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II.  MICROPULSATIONS  AND  THE  IONOSPHERE 


A.  THE  PREDICTION  PROBLEM 

Of  fundamental  concern  is  the  size  of  the  effect  being  studied. 

How  great  a  fluctuation  in  the  ionospheric  electron  concentration  is  to 
be  expected,  and  how  much  influence  will  such  a  fluctuation  have  on  a 
probing  HF  wave? 

Starting  with  the  theory  of  Pc  1  propagation  outlined  in  Chapter  I 
and  making  simplifying  assumptions  as  necessary,  we  can  proceed  to  make 
approximate  calculations  of  the  magnitudes  in  question.  The  investi¬ 
gation  will  be  concerned  with  the  middle-latitude  nighttime  ionosphere 
in  the  northern  hemisphere,  and  the  presence  of  Pc  1  micropulsation 
energy  travelling  southward  will  be  assumed.  The  approach  will  be  to 
use  the  simplest  models  that  account  for  the  relevant  physical  proper¬ 
ties  of  this  system. 

In  what  follows,  whenever  the  coordinate  system  needs  to  be  related 
to  the  earth,  the  directions  x,  y,  z  refer  to  south,  east  and  up 
respectively.  The  geomagnetic  field  is  assumed  to  be  oriented  parallel 
to  the  z  axis  and  thus  to  be  vertical.  Note  that  the  discussion 
initially  refers  to  a  quite  general  plasma  ;  parameters  are  then  chosen 
so  that  this  plasma  models  the  physical  earth- ionosphere  system. 

B.  IONOSPHERIC  ELECTRON-DENSITY  PERTURBATION 

For  the  purpose  of  estimating  th«  size  of  the  ionospheric  electron- 
density  perturbation  we  will  examine  the  simplest  plasma  that  will  , 
permit  propagation  of  a  hydromagnetic  wave.  Consider  an  unbounded, 
cold  and  collisionless  plasma  which  in  the  unperturbed  state  is  homo¬ 
geneous  and  at  rest.  We  require  that  this  plasma  be  pervaded  by  a 

uniform  magnetostatic  field  of  strength  B  tesla,  and  then  establish  a 

o 

cartesian  coordinate  system  in  the  plasma  so  oriented  that  the  z  axis 
is  parallel  to  the  magnetostatic  field.  Then  all  of  the  possible  modes 
of  propagation  of  any  monochromatic  plane  electromagnetic  wave  with  its 
wavr  normal  lying  in  the  x-z  plane  are  described  by  the  equation 
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2  2  2 

Si  -n  cos  8  -ie  n  sin9cos9 

E 

x  X 

X 

2 

ie  Si-n  0 

E 

X 

y 

2  2  2 
n  sinfjcosQ  0  ®ll  “n  sin  9 

E 

.  z 

(1) 


Here  n  is  the  refractive  index  of  the  plasma  and  0  is  the  angle  between 

the  wave  normal  and  the  positive  z-axis.  The  vector  E  =  (E  E  E  ) 

x  y’  z 

represents  the  (time-varying)  electric  field  of  the  electromagnetic  wc'e 
and  the  e's  are  components  of  the  effective  tensor  permittivity  of  the 
plasma  as  given  in  Appendix  C,  where  equation  (1)  is  derived. 

Physically  we  are  interested  in  a  hydromagnetic  wave  travelling 
horizontally  in  the  mid-latitude  ionosphere  and  therefore  crossing  lines 
of  the  earth's  magnetostatic  field  at  steep  incidence.  We  will  simulate 
this  situation  in  our  model  by  considering  a  wave  moving  along  the 
positive  x-axis.  When  (1)  is  solved  with  6  set  equal  to  90°,  it  can  be 
shown  (Clemmow  and  Dougherty,  1969)  that  there  are  two  possible  modes 
in  which  such  a  wave  can  propagate.  However,  in  a  plasma  having  an 
electron  density  approximating  that  of  thJ  ionospheric  F  region  and  at 
the  low  frequencies  of  interest  in  this  study  (that  is,  frequencies  near 
1  Hz),  only  one  of  these  modes  can  propagate.  For  this  mode 


Si 


When  this  value  of  n  is  inserted,  (1)  reduces  to 

e 

Ex  =  1  Ey  and  Ez  =  0  • 

Invoking  Maxwell's  equation  for  the  curl  of  E,  the  components  of  the  wave 
magnetic  field  B  are  discovered  to  be 

B  =  -  E  and  B  =  B  =  0  , 
z  c  y  x  y  ' 

where  c  is  the  vacuum  speed  of  light. 
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Since  a  wave  frequency  of  1  Hz  is  much  less  than  the  ion  cyclotron 
frequency  of  our  ionosphere-approximating  plasma,  the  following  approx¬ 
imations  nry  be  made  as  discussed  in  Appendix  C: 

2 

c 

e.fcs  —  and  e  0 

2  x  ’ 

VA 

where  va  is  the  Alfven  speed.  The  wave  field  components  now  reduce  to 

Ey  =  Vz  and  EX  =  EZ  -  BX  =  By  "  0  * 

We  now  confine  our  attention  to  the  electron  motions  controlled  by 
these  fields.  Let  N  and  AN  be  respectively  the  unperturbed  electron 
number  density  and  the  perturbation  thereof  caused  by  the  passage  of  our 
wave.  Then  from  the  charge  conservation  equation, 


AN 


where  Vx  is  the  x  component  of  electron  velocity.  From  the  Lorentz 
force  equation 


These  last  two  equations  combine  to  give 


AN  \ 

N  ~  B  (2) 

o 


That  is,  the  fractional  electron-density  perturbation  is  equal  to  the 
ratio  of  the  amplitude  of  the  wave  magnetic  field  component  parallel  to 
the  magnetostatic  field  to  the  amplitude  of  the  magnetostatic  field 
itself.  This  equation  is  the  central  relationship  for  this  experiment 
since  it  enables  calculation  of  the  electron-density  perturbation 
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when  the  wp.ve  magnetic  field  strength  is  known,  and  vice  versa, 

C.  FIELD  VARIATION  WITH  ALTITUDE 

The  plasn.i  model  discussed  in  Section  B  is  unbounded  and  uniform  - 
that  is,  it  exists  in  some  imaginary  space  and  has  no  connection  with 
the  earth.  This  is  a  reasonable  model  for  the  ionosphere  but  cannot, 
however,  give  much  idea  of  what  happens  at  the  earth's  surface.  And  a 
connection  with  the  earth  must  be  made  since  the  only  method  we  have  at 
present  of  estimating  the  size  of  micropulsation  wave  fields  in  the 
ionosphere  is  by  making  measurements  at  the  ground. 

There  are  two  principal  objections  to  the  model  used  in  Section  B. 
The  first  is  that  the  actual  electron  density  in  the  atmosphere  is  not 
uniform  with  altitude,  bat  varies  in  a  manner  similar  to  that  shown  in 
Figure  10  of  Appendix  B.  A  height-dependent  electron  density  must  be 
introduced  if  we  are  to  modol  waves  at  the  ground  and  in  the  ionosphere 
simultaneously.  Secondly,  there  is  actually  a  sharp  boundary  at  the 
surface  of  the  earth  since  the  earth  is  a  good  electrical  conductor  at 
frequencies  near  1  Hz.  Our  model  should  be  similarly  bounded  if  effects 
near  the  existing  real  boundary  are  to  be  approximated. 

A  convenient  model  which  meets  the  requirements  laid  down  in  the 
previous  paragraph  is  one  used  by  Greifinger  and  Greifinger  (1968) .  A 
plot  of  the  electron-density  variation  with  height  for  this  model  is 
shown  in  Figure  2.  The  coordinate  system  is  related  to  the  earth  in  the 
manner  discussed  in  Section  A.  The  earth  itself  is  taken  to  be  flat  and 
perfectly  conducting,  providing  a  sharp  bound  at  z  =  0.  Hence  the  model 
is  semi-infinite  in  the  z  direction  but  is  still  unbounded  in  the  x  and 
y  directions.  The  dimensions  of  the  model  are  consistent  with  nighttime, 
sunspot-minimum  conditions. 

The  Greifingers  develop  equations  in  their  article  which  describe 
the  fields  associated  with  a  1-Hz  hydromagnet ic  wave  travelling  horizon¬ 
tally  in  this  model.  The  magnetic  field  components  of  such  u  wave 
moving  southward  are  plotted  in  Figure  2.  Now  Bz  in  the  ionosphere  (Bzm) 
is  proportional  to  the  fractional  electron-density  perturbation  there  to 
the  extent  that  equation  (2)  can  still  be  considered  valid;  moreover, 

B  at  the  ground  (B  )  can  be  measured.  Since  the  ratio  of  these  two 
x  xo 
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HEIGHT  (km) 


I iguie 


ELECTRON  DENSITY  (electrons/m3) 


2. 


ion- Phere 

tW°  COmPonen^  ot  the  magnetic  wave  ulld'TT 
micropulsation  travelling  in  this  ionosphere  B  at  tt 
ground  is  arbitrarily  chosen  to  be  100  mV  The  rati 
B  in  the  ionosphere  tn  r  q*  ' '  stio  o 

z  nospnere  to  Bx  at  the  ground  is  roughly  2-5 
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quantities  is  a  constant  (for  a  fixed  wave  frequency),  knowing  its  value 

enables  us  tc  estimate  the  size  of  the  wave  fields  in  the  ionosphere,  and 

thus  of  the  electron-density  perturbation  there,  given  a  measurement 

made  on  the  ground.  Examination  of  Figure  2  reveals  that  B  /B  has  a 

zm  xo 

value  of  roughly  2.5  for  our  model  at  1  Hz. 


D.  DFTECTING  THE  DISTURBANCE 

Given  that  there  is  a  certain  electron-density  disturbance  in  the 
ionosphere,  how  is  the  ph^se  of  a  probing  HF  signal  affected?  To  answer 
this  question  let  us  again  consider  the  model  discussed  in  Section  C, 
except  that  we  now  set  B0  =  0  since  the  presence  of  the  magnetostatic  field 
does  not  affect  the  essence  of  the  HF  propagation  now  under  discussion. 

If  N(z)  is  the  electron-density  distribution  with  height,  then  the 
square  of  the  refractive  index  is  given  by 


1  /  2 

Here  f  (=  9N  '  in  mks  units)  is  the  electron  plasma  frequency,  f  is 
P  2 

the  HF  wave  frequency,  and  A  =  8l/f  is  a  convenient  parameter  (see 
Budden,  1961)  .  Using  this  equation  we  can  solve  for  the  refractive 
index  perturbation  An  caused  by  an  electron  density  change  AN  : 


2  !/2 
An  =  -n  +  (n  -A AN) 


(3) 


The  next  step  is  to  examine  the  phase  path  of  a  vertically-incident 

radio  wave.  If  z  is  the  height  at  which  the  wave  is  reflected,  then  the 
r 

phase  path  (see  Budden)  is 


P 


z 

r 


2  / 
0 


n  dz 
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For  refractive  index  perturbation  An,  and  assuming  that  the  resulting 
Az^  is  small,  the  phase-path  change  is 

z 

r 

AP  »  2  /  An  dz  (4) 

0 

When  the  phase-path  change  and  the  free-space  HF  wavelength  are  known, 
then  the  number  of  degrees  of  phase  change  of  the  reflected  radio 
wave  is  easily  calculated. 

E.  A  NUMERICAL  ESTIMATE 

To  get  a  quick  numerical  estimate  of  the  phase  change  that  might 

result  from  the  passage  of  a  typical  Pc  1  micropulsation  in  the 

ionosphere,  let  us  evaluate  equation  (4)  using  the  model  of  Figure  2. 

2  4  2 

Choose  =  400  km  and  f  =  —  f  .  Assuming  that  equation  (2)  holds 

everywhere  in  the  ionized  region,  and  using  the  result  of  Section  C 

that  B  /B  rs  2,5  at  1  Hz, 
zm  xo 

AN  2.5  x  100  mV  _  -6 

—  »  - - -  =  5  X  10 

N  0.5  gauss 


for  a  measured  field  on  the  ground  of  100  mV  peak-to-peak  and  a  static 
field  of  0.5  gauss.  We  will  take  AN/N  to  be  constant  in  the  region  of 
interest.  Substituting  into  (3)  and  (4), 

400  km 

AP  =  -7.5  x  10  /  dz  =  -0.75  m 

300  km 

For  an  HF  wavelength  of  72  meters,  this  represents  0.75  x  360° /72  «  3.8° 
of  phase  sung.  In  other  words,  we  might  reasonably  expect  that  the 
phase  of  the  probing  signal  would  vary  over  a  range  of  something  like  a 
few  degrees  in  a  typical  situation. 
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III.  NOISE 


A.  NOISE  SOURCES 

The  problem  of  detecting  a  small  phase  modulation  of  an 
ionospherically-ref lected  wave  is  complicated  by  the  presence  of  a 

number  of  contaminating  signals.  These  signals  can  be  roughly  classified 
in  the  following  manner  : 

(i)  background  noise, 

(ii)  phase  noise  imposed  by  the  ionosphere, 

(iii)  spurious  reflections  of  the  transmitted  signal. 

Interference  from  other  transmitters  is  neglected  since  the  receiving 
system  is  narrowband  and  unwanted  carriers  are  easily  avoided.  Noise 
introduced  by  the  detection  system  is  generally  negligible  compared  to 
the  external  noise  and  is  not  discussed  here. 

B,  BACKGROUND  NOISE 

The  background  noise  is  assumed,  not  unreasonably,  to  be  Gaussian 
white  noise.  This  white  noise  is  filtered  by  the  bandpass  characteristic 
of  the  receiver  and  becomes  narrowband  noise  which  can  be  represented  by 

the  expression  n(t)=an  cos  (wt+cpn)  where  an  and  cpn  are  random  processes 
(Carlson,  1968) 

The  noise  adds  to  the  desired  signal  x(t)  =  a  cos  (wt  +  cp> 
to  yield 

s(t)  =  x ( t )  +  n(t) 

=  a  cos  (ut  +  <p)  +  a  cos  (wt  +  cc  ) 

n  Tn 

=  M  cos  (tot  +  ¥) 

where  expressions  f>r  M  and  ¥  can  be  found  by  algebraic  manipulation. 

When  the  assumption  is  made  that  an  «  a  most  of  the  time,  as  is 
usually  the  case  in  practice,  then  ¥  is  given  by 

a 

¥  =  cp  -  sin  (cp  -  cpn)  . 


Preceding  page  blank 
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The  power  spectrum  of  the  second  term  in  this  expression  is  just  the 
power  spectrum  of  th'3  original  noise  but  diminished  by  the  factor  l/a2 
(Carlson,  1968).  Hence  the  additive  noise  is  transformed  directly  to 
additive  Gaussian  narrowband  phase  noise,  but  with  the  noise  power 
reduced  according  to  the  strength  of  the  transmitted  carrier.  For  the 
experiment  under  consideration  in  this  report,  the  carrier  and  noise 
levels  were  generally  such  that  background  noise  was  not  a  limiting 
factor . 

C.  IONOSPHERIC  PHASE  NOiSE 

The  electron  density  and  thus  the  refractive  index  along  the  path 
of  a  radio  wave  travelling  through  the  ionosphere  are  continually  and 
randomly  varying  with  time.  This  variation  modulates  the  phase  of  our 
received  wave  and  hence  is  a  noise  source.  The  effect  is  probably  quite 
small,  however,  especially  since  the  expression  for  the  received  phase 
involves  an  integral  over  the  path.  In  any  event,  such  noise  cannot 
be  easily  separated  from  that  due  to  spurious  reflections  and  will  be 
neglected  from  here  on. 

D.  SPURIOUS  REFLECTIONS 

Energy  radiated  from  a  simple  horizontal  dipole  antenna  can  reach 
a  receiving  antenna  some  short  distance  away  by  a  multitude  of  paths. 
Assuming  that  the  groundwave  is  blocked  in  some  way,  all  the  signal  paths 
involve  ionospheric  reflection.  The  number  of  paths  and  their  propaga¬ 
tion  characteristics  depend  on  the  roughness  and  turbulence  of  the 
ionosphere  at  any  particular  time. 

The  simplest  case  that  one  can  expect  to  encounter  occurs  at  times 
when  the  ionosphere  is  effectively  a  smooth  reflector  of  radio  waves. 

In  such  a  situation  the  signal  arriving  at  the  receiver  is  the  combination 
of  several  components  since  the  original  wave  has  divided  into  ordinary 
and  extraordinary  components  and  because  some  of  the  wave  energy  has  made 
multiple  ground-ionosphere-ground  hops  before  being  received.  The 
ordinary  and  extraordinary  components  reflect  at  ditferent  heights  but 
have  comparable  amplitudes  and,  assuming  the  presence  of  a  periodic 


SEL-72-021 


14 


electron-density  perturbation,  comparable  levels  of  phase  modulation 
Each  successive  hop,  however,  has  a  greater  level  of  phase  modulation 
than  itr  iredecessor  because  it  has  passed  through  the  modulating  region 
more  often,  but  is  lower  in  amplitude  because  it  has  suffered  more 
ionospheric  attenuation  and  has  travelled  farther.  Each  component  is 
in  general  reflected  from  a  different  point  in  the  ionosphere,  but  these 
points  are  within  a  few  tens  of  kilometers  of  each  other  for  near- 
vertical  incidence.  Since  we  are  interested  in  ionospheric  distuioances 
with  large  horizontal  scale  size  (on  the  order  of  1000  km),  we  can 
assume  that  the  modulation  imposed  on  each  component  is  approximately 
in  phase  with  the  modulation  imposed  on  all  the  other  components. 

However,  the  carrier  phase  of  each  component  is  completely  random  with 
respect  to  the  carrier  phases  of  the  other  components,  in  addition, 
the  carrier  frequency  of  each  component  has  a  distinct  doppler  shift 

since  each  reflection  point  has  in  general  a  distinct  effective  vertical 
velocity. 

Some  "./uputer  simulations  were  made  of  the  smooth-ionosphere 
situation  and  the  results  of  these  indicate  that  doppler  shifts  of  the 
component  carrier  frequencies  provide  most  of  the  noise  in  the  phase 
spectrum  of  the  received  signal.  The  net  phase  modulation  level  at 
the  receiving  antenna  is  not  very  different  from  the  level  that  would 
be  expected  if  there  were  only  one  path  between  transmitter  and  receiver, 
but  the  phase  spectrum  fills  with  noise  as  doppler  shifts  are  introduced! 
However,  for  moderate  doppler  shifts  consistent  with  reasonably  settled 
ionospheric  conditions,  phase  modulation  levels  of  a  few  degrees  of  phase 
swing  can  be  readily  detected  above  the  noise;  Of  course  when  the 
ionosphere  becomes  rough  and  turbulent,  then  the  number  of  signal 
components  and  the  size  of  the  doppler  offsets  increase  and  more  phase 
noise  is  produced,  eventually  obscuring  the  phase  modulation.  In 
practice,  there  is  some  "normal”  level  of  ionospheric  disturbance  which 
determines  a  minimum  level  of  phase  modulation  that  can  be  detected  ; 
this  threshold  is  determined  by  experimental  measurement. 
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IV.  THE  APPARATUS 


A.  EXPERIMENTAL  ARRANGEMENT 

The  experiment  was  implemented  using  a  simple  continuous-wave  (CW) 
transmitter  in  conjunction  with  a  narrowband  receiver  and  phase 
demodulator.  The  transmitter  site  was  established  at  a  location  on  the 
Pacific  coast  of  California  about  28  km  irom  the  receiver  site  on  the 
Stanford  University  campus.  The  geographical  layout  is  illustrated  in 
Figure  3.  The  range  of  h  ,11s  running  between  the  two  sites  eliminated 
groundwave  propagation.  Hence  the  main  signal  path  between  the  two 
sites  was  the  one-hop  ionospheric  route.  This  was  essentially  a  vertical 
path  due  to  the  relatively  small  separation  between  transmitter  and 
receiver . 

The  transmitter  radiated  approximately  100  watts  at  about  4  MHz 
from  a  horizontal  half-wave  dipole  antenna.  A  similar  antenna  was  used 
at  the  receiver  site.  Call-letter  keying  and  a  calibration  modulation 
were  applied  to  the  transmitted  signal  as  described  in  Section  C  of  this 
chapter , 

The  receiving  system  is  outlined  in  Figure  4.  The  4-MHz  input  is 
down-converted  to  100  Hz  with  a  bandwidth  of  10  Hz.  From  this  point 
the  signal  is  demodulated,  recorded,  and  eventually  processed  as  indi¬ 
cated  by  the  succeeding  blocks  in  Figure  4. 

Supplementing  this  primary  equipment,  a  geomagnetic  micropulsation 
recording  station  and  a  vertical  ionosonde  were  in  operation  on  the 
Stanford  campus.  The  micropulsation  station  continuously  recorded  the 
time-varying  horizontal  north-south  (geographic)  component  of  the 
earth's  magnetic  field.  The  ionosonde  provided  vertical  ionograms 
several  times  each  hour  of  the  day. 


B .  THE  DEMODULATOR 


The  demodulator  is  a  phase-locked  loop  as  diagrammed  in  Figure  5. 
The  voltage  controlled  oscillator  (VC0)  operates  nominally  at  the  same 
frequency  and  phase  as  the  input  signal.  It  follows  slow  input-phase 
variations  and  maintains  the  phase— detector  outputs  at  zero.  However, 
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Figure  3.  Geographical  layout  of  the  experimental  apparatus 
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Figure  5.  Block  diagram  of  the  demodulator. 


relative  amplitude 


Figure  6.  Frequency  response  of  the  receiving  system  to  phase  modu¬ 
lation  . 
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due  to  the  effect  of  the  lo.pass  amplifier,  the  VCO  cannot  fotlo.  rapid 

input-phase  variations.  Hence  rapid  variations  .ill  appear  as  output  of 

both  phase  detectors.  In  other  words,  the  loop  Itself  has  a  hlghpass 

response  to  phase  modulation.  However,  the  narrow  bandwldths  of  the 

preceding  stages  of  the  receiver  limit  the  high-frequency  response  of  the 

system  and  the  overall  effect  is  to  produce  a  bandpass  response  to  phase 

modulation.  The  measured  response  of  the  receiver  system  is  shown  in 
Figure  6. 


C.  CALIBRATION  SIGNAL 

Noise  is  a  major  component  of  the  demodulator  output.  There  are 
contributions  from  both  the  ionosphere  and  the  receiver.  Since  the 
signals  being  sought  are  very  small,  it  is  necessary  to  know  at  any  time 
•hat  the  smallest  detectable  signal  is.  For  this  reason  a  calibrated 
Phase  modulation  is  applied  to  the  transmitted  signal.  intensity  is 

stepped  over  a  range  of  values  in  order  that  the  minimum  detectable 
signal  at  any  point  in  time  can  be  determined  simply  by  examining  the 
output  record.  The  standard  calibration  format  is  given  in  Table  ,. 

The  receiving  system  is  capable  of  resolving  less  than  0.5-  of  phase 
swing,  but  this  level  of  performance  is  no,  attainable  under  actual 
operating  conditions  because  of  the  presence  of  ionospheric  noise 
(as  described  in  Chapter  III). 

D.  RECORDING  AND  REDUCTION 

The  demodulator  output  is  recorded  on  magnetic  tape  and  on  a  strip 
chart.  The  chart  records,  examples  of  which  are  shown  in  F,gure  7, 
give  a  real-time  determination  of  propagation  conditions  and  general 
system  operation.  The  magnetic  tape  provides  the  main  data  record 

It  is  often  very  difficult  to  understand  amplitude-time  records 
because  of  the  complexity  of  the  signals  and  the  presence  of  noise. 

Hence  the  data  stored  on  magnetic  tape  is  converted  to  frequency-time 
form  for  easier  study.  This  is  accomplished  using  an  analog  spectrum 
analyzer  and  a  tape  playback  speed  thirty-two  times  greater  than  the 
recording  speed.  The  increased  speed  keeps  the  amount  of  time  required 
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TABLE  1 


Modulation 

and  Keying  Schedule 

Minutes  after 

Hour  or  Half  Hour 

Modulation  Level 
(peak-to-paak  phase  swing) 

0 

call-letter  keying 

0.0  - 

1.5 

carrier  off 

1.5  - 

3.0 

0° 

3.0  - 

7.5 

32° 

7.5  - 

12.0 

16° 

12.0  - 

16.5 

8° 

16.5  - 

21.0 

4° 

21.0  - 

25.5 

2° 

25.5  - 

28.5 

0° 

28.5  - 

30.0 

carrier  off 

30 

call-letter  keying 
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required  for  playback  at  a  practical  level  (about  eleven  minutes  for 
six  hours  of  data)  and  provides  a  signal  frequency  more  suited  to  the 
spectrum  analyzer. 

The  analyzer  output  is  recorded  either  photographically  or  by  a 
facsimile  process.  In  each  case  the  record  shows  a  plot  of  frequency 
versus  time  of  day,  with  signal  strength  indicated  by  relative  darkening 
of  the  film  negative  or  facsimile  paper.  The  film  record  in  produced  by 
pulling  film  across  the  face  of  an  oscilloscope  whose  x  and  z  axes  are 
driven  respectively  by  the  frequency  and  amplitude  outputs  of  the 
analyzer.  The  facsimile  record  is  formed  on  chemically-treated  paper 
by  a  special  recorder.  An  advantage  of  the  facsimile  process  is  that 
it  provides  immediately-visible  output  which  allows  for  the  expeditious 
setting  of  gains  and  speeds,  whereas  film  must  first  be  developed  before 
any  feedback  can  be  applied.  An  example  of  the  facsimile  output  is 
shown  in  Figure  8. 
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V .  RESULTS 


A.  EXPERIMENTAL  PROCEDURE 

The  nominal  procedure  was  to  record  phasa  variations  of  the  HF  signal 
for  a  six-hour  period  each  night  (usually  from  2300  to  0500  PST).  This 
was  done  intermittently  over  the  interval  from  .June  18,  1970  to 
December  31,  1970,  and  essentially  without  interruption  from  January  1, 
1971  to  June  25,  7.971.  Equipment  problens  were  responsible  for  most 
nights  missed.  The  six-hour  recording  time  was  dictated  by  properties 
of  the  tape  recorder,  but  proved  long  enough  to  provide  coverage  of  most 
Pc  1  micropulsation  events. 

Any  time  that  the  micropulsation  records  showed  that  a  Pc  1  event 
had  occurred  on  the  ground,  the  phase  records  were  searched  for  evidence 
of  coincident  fluctuations.  The  ionograms  for  such  times  were  also 
consulted  to  determine  the  proximity  of  the  operating  frequency  to  the 
ionospheric  critical  frequency. 

B.  Pc  1  MICROPULSATION  OCCURRENCES 

Over  the  year  that  the  experiment  was  in  operation,  Pc  1  micro¬ 
pulsation  activity  occurred  (at  the  ground)  on  an  average  of  5.5  nights 
per  month.  This  activity  is  plotted  on  a  monthly  basis  in  Figure  9. 

A  noticeable  feature  of  the  plot  is  the  fall-winter  minimum  and  spring- 
summer  maximum  of  occurrences  . 

C.  EXPERIMENTAL  RESULTS 

During  the  course  of  the  experiment  there  were  only  9  nights  on 
which  reasonably-sized  micropulsations  (50  mV  or  greater  peak-to-peak 
amplitude)  occurred  while  the  phase-measuring  equipment  was  operating 
properly.  No  coincident  phase-path  fluctuations  could  be  detected  at 
any  of  these  times. 

Some  of  these  events,  including  the  largest  ones,  are  listed  in 
Table  2  along  with  values  of  relevant  parameters.  The  probe  frequency 
in  all  cases  was  about  4  MHz.  The  detection  ratio  used  in  Table  2  is 
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Some  Experimental  Numbers 
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the  ratio  of  the  measured  micropulsation  amplitude  to  the  minimum 
detectable  phase-modulation  level,  and  is  a  measure  of  the  likelihood 
of  observing  phase-path  fluctuations.  Larger  micropulsations  or  a 
lower  detection  threshold  give  higher  values  for  this  ratio. 

In  Chapter  II  the  ratio  Bz  /Bxq  was  estimated  to  be  about  2.5 
tor  1-Hz  micropulsations.  Noting  that  the  largest  micropulsations 
observed  had  a  peak-to-peak  amplitude  on  the  ground  (B^q)  200  mV* 

we  can  estimate  that  the  fractional  electron-density  perturbation  in 
the  ionosphere  must  have  been  less  than  about  10  (taking  An/N=Bz.  Bq 

and  B  =0.5  gaust) . 

o 

Table  2  demonstrates  that  the  normal  noise  level  corresponds  to 
a  minimum  detectable  phase  swing  of  about  4° .  Since  the  background 
noise  level  was  low  enough  to  be  neglected,  this  noise  presumably 
results  from  the  effects  of  spurious  reflections  as  was  discussed  in 
Chapter  III. 

D.  DISCUSSION 

The  phase  swing  estimated  for  a  typical  micropulsation  in  Section  E 
of  Chapter  II  is  of  the  s..me  order  of  magnitude  as  the  minimum 
detectable  phase  swing  measured  during  the  experiment.  Although  this 
leads  to  the  expectation  that  the  largest  micropulsations  might  have 
produced  an  observable  disturbance,  it  is  not  surprising,  in  light  of 
the  various  approximations  made  in  the  preceding  pages,  that  no  such 
disturbances  were  indeed  detected. 

The  models  used  in  Chapter  II  neglect  the  inclination  of  the 
magnetostatic  field  and  do  not  provide  a  very  realistic  N(z)  profile. 

The  assumption  is  made  that  equation  (2)  is  valid  at  all  heights,  and 
only  1-Hz  micropulsations  are  considered.  The  effect  of  the  bandpass 
characteristic  of  the  phase-locked  loop  is  ignored,  as  is  the  inte¬ 
gration  time  of  the  spectrum  analyzer.  All  of  these  matters,  and 
probably  others  as  well,  would  have  to  be  taken  into  consideration  if 

a  complete  description  of  the  processes  involved  in  this  experiment 
* 

were  to  be  given . 

Much  of  this  work  has  since  been  done  by  the  author,  but  is  not 
included  here  because  it  does  not  materially  affect  the  results  of  this 
study . 
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The  failure  of  the  experiment  to  detect  micropulsation  disturbances 
reduces  to  the  statement  that  the  signal  level  must  have  been  less 
than  the  detection  threshold  of  the  receiving  system.  The  sensitivity 
of  the  equipment  might  be  improved  somewhat  by  such  measures  as  using 
polarized  and/or  directive  antennas  to  cut  down  on  the  spurious  reflec¬ 
tions  received,  or  by  improving  the  time-  and  frequency-response 
characteristics  of  the  detection  and  data-reduction  systems.  However, 
assuming  that  little  improvement  can  be  made  in  the  detection  scheme, 
the  best  hope  of  achieving  some  positive  results  might  be  to  relocate 
the  experiment  at  a  more  northerly  latitude  where  the  amplitude  of  the 
micropulsations  can  be  expected  to  be  larger. 
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VI .  CONCLUSIONS 


Based  on  the  approximate  calculations  made  in  Chapter  II,  and  using 
the  amplitude  of  the  largest  micropulsations  observed  on  the  ground,  it 
can  be  estimated  that  the  fractional  electron-density  perturbation  pro¬ 
duced  in  the  ionosphere  by  travelling  Pc  1  micropulsations  is  normally 

less  than  about  10  parts  per  million  for  geomagnetic  latitudes  lower  than 
45°  . 

In  addition,  the  experimental  detection  scheme  described  in  this 
report  is  limited,  probably  by  the  effect  of  spurious  reflections,  to 
a  minimum  detectable  phase  modulation  level  of  about  4°  of  phase  swing 
at  night. 

Working  within  the  constraints  of  an  upper  bound  on  the  electron- 
density  perturbation  and  a  minimum  detectable  phase  change,  the 
experiment  was  unable  to  detect  any  ionospheric  disturbances  associated 
with  Pc  1  micropulsations. 
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APPENDIX  A 


MICROPULSATION  CLASSIFICATION  AND  UNITS 
1.  CLASSIFICATION 

Based  on  experimental  findings,  geomagnetic  micropulsations  are 
divided  into  two  main  classes  -  those  of  a  regular  and  mainly  continuous 
character,  and  those  with  an  irregular  pattern.  The  first  class  covers 
the  whole  range  of  micropulsations  with  periods  from  about  0.2  to  600 
seconds,  and  is  divided  into  subgroups  based  on  physical  and  morpho¬ 
logical  properties.  Table  3  shows  the  classification  and  notation  for 
these  continoes  micropulsations. 

Table  3 

Continuous  Micropulsations 

Period  Range 


Notation 

(seconds) 

Pc  1 

0.2  -  5 

Pc  2 

5-10 

Pc  3 

10  -  45 

Pc  4 

45  -  150 

Pc  5 

150  -  600 

The  second  main  class  of  micropulsations  is  characterized  by  their 
irregular  form,  their  close  connection  with  disturbances  of  the  geo¬ 
magnetic  field,  and  their  correlation  with  various  upper  atmospheric 
phenomena.  For  this  class  there  are  two  subgroups  as  shown  in  Table  4. 

Table  4 

Irregular  Micropulsations 

Period  Range 

Notation  I  seconds! 

Pi  1  1-40 

Pi  2  40  -  150 
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This  description  is  taken  from  Jacobs  (1970)  where  further  details 
may  be  found. 

2.  UNITS 

The  earliest  workers  in  the  field  of  magnetism  envisioned  the  magnetic 
field  in  terms  of  lines  of  force  emanating  from  magnetic  poles,  where  a 
pole  was  considered  to  be  the  magnetic  equivalent  of  an  electric  point 
charge.  Although  magnetic  poles  have  never  been  found  to  exist  as 
independent  entities,  the  concept  of  lines  of  force  remains  useful,  and 
in  fact  the  magnetic  field  strength  B  is  conveniently  described  in  terms 
of  such  lines.  (The  force  that  we  are  talking  about  here  is  that  which 
would  be  exerted  on  a  small,  positive  magnetic  "charge"  placed  in  the 
field. )  At  any  point,  B  is  directed  along  the  force  lines  end  is  pro¬ 
portional  in  magnitude  to  the  line  density  (that  is,  the  number  of  lines 
passing  through  unit  cross-sect±onal  area).  Note  that  B  is  histori¬ 
cally  called  the  magnetic  induction,  but  is  termed  the  magnetic  field 
strength  in  this  report  since  this  is  felt  to  be  a  more  descriptive  name. 

The  unit  of  B  in  the  mks  system  is  the  tesla  (or  weber  /m2)  .  For 
convenience,  and  following  historical  precedent,  smaller  units  are 
usually  employed,  namely  the  gauss,  the  gamma,  and  the  mtiligamma. 

Table  5  shows  the  relationships  among  these  units. 

Table  5 

Units  of  Magnetic  Field  Strength 
-4  2 

1  gauss  (G)  =  10  weber/m  =  100  microteslas 

1  gamma  ( y)  =  10  5  gauss  =  l  nanotesla 

1  milligamma  (mY)  =  10~8  gauss  =  1  picotesla 
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APPENDIX  B 


DEFINITIONS 

The  following  information  is  taken  in  large  measure,  and  often 

verbatim,  from  Clemmow  and  Dougherty  (1969),  Rishbeth  and  Garriott 
(1969),  and  Utlaut  and  Cohen  (1971)  . 

1 .  PLASMA 

A  Plasma  is  any  assemblage  of  particles,  some  or  all  of  which  are 
charged,  that  can  from  the  macroscopic  view  be  regarded  as  a  gas  which 
in  its  equilibrium  state  is  everywhere  locally  neutral.  The  fact  that 
individual  particles  of  the  gas  are  charged,  and  therefore  generate  and 
interact  with  electromagnetic  fields,  gives  rise  to  many  distinctive 
phenomena . 

The  plasma  frequency  is  the  frequency  of  natural  oscillation 
associated  with  the  electrostatic  restoring  force  arising  from  a  small 
displacement  of  charge  in  a  plasma.  For  plasmas  in  a  magnetic  field, 
it  is  often  necessary  to  consider  the  gyro  frequency  (or  cyclotron 
^frequency)  ;  this  is  the  frequency  at  which  a  charged  particle  moving  at 
right  anglos  to  the  field  is  forced  to  gyrate  about  the  field  lines. 

When  relative  bulk  motion  of  two  plasmas  is  considered,  the  possi¬ 
bility  of  plasma  instability  arises.  That  is,  electromagnetic  waves 
can  be  generated  or  amplified  in  one  of  the  plasmas  by  the  process  of 

extracting  energy  from  the  relatively-moving  particles  of  the  other 
plasma . 

2.  HYDROMAGNET I C  WAVE 

^ow-frequency  electromagnetic  waves  (that  is,  waves  with 
frequencies  much  less  than  the  ion  gyro  frequency)  travelling  through  a 
plasma  in  a  magnetostatic  field  have  the  property  that,  considering 
the  associated  motion  of  the  charged  particles,  the  velocity  component 
perpendicular  to  the  magnetostatic  field  is  approximately  the  same  for 
the  electrons  as  for  the  ions.  In  its  motion  across  the  magnetostatic 
field,  the  plasma  therefore  moves  as  a  whole.  Furthermore,  these  waves 
travel  at  a  speed  which  depends  only  on  the  mass  density  of  the  plasma 
and  the  strength  of  the  magnetostatic  field;  this  is  the  Alfven  speed. 
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Such  waves  are  variously  called  hydromagnetic ,  magnetohydrodynamic  (mhd) 
or  Alfven  waves. 

3 .  IONOSPHERE 

The  ionosphere  is  that  region  of  the  earth's  atmosphere  where  ions 
and  electrons  are  present  in  sufficient  quantities  to  affect  the  propa¬ 
gation  of  radio  waves.  This  region  begins  at  an  altitude  of  about  50  km 
and  extends  upward  as  shown  by  plots  of  electron  concentration  versus 
altitude  such  >.s  the  one  given  in  Figure  10.  The  shape  of  these  profiles 
varies  diurnally,  seasonally,  with  sunspot  number,  and  with  latitude  and 
longitude.  The  profile  peak  usually  occurs  between  200  and  400  km;  the 
portion  of  the  ionosphere  in  the  vicinity  of  the  peak  is  commonly  termed 
the  F  region  or  layer.  In  the  F  region  the  peak  densities  (number  densi¬ 
ties  or  concentrations)  of  ions  (predominantly  0+  ions)  and  electrons  are 
each  about  1012  particles  per  cubic  meter.  The  neutral  particle  density 
is  about  1015  per  cubic  meter,  implying  only  0.1  percent  or  less  ioniza¬ 
tion  in  this  region.  At  about  1500  km,  however,  the  ion  and  neutral- 
particle  densities  become  roughly  comparable,  and  this  is  a  convenient 
point  to  label  as  the  upper  boundary  of  the  ionosphere. 

The  free  electrons,  rather  than  the  ions,  interact  more  strongly 
with  radio  waves  since  they  are  much  less  massive,  and  it  is  the  electron 
plasma  frequency  of  the  ionosphere  which  determines  the  level  at  which 
impinging  radio  waves  are  reflected.  The  maximum  plasma  frequency  (which 
occurs  at  the  electron  density  peak)  is  called  the  critical  or  penetration 
frequency  of  the  ionosphere  since  vertically-incident  waves  of  higher 
frequency  will  not  be  reflected  but  will  pass  through  the  ionosphere. 

The  critical  frequency  typically  lies  between  3  and  6  MHz  at  night,  and 
between  6  and  12  MHz  during  the  day. 

4 .  MAGNETOSPHERE 

The  magnetosphere  is  that  region  of  the  earth’s  atmosphere  where 
the  particle  dynamics  are  primarily  controlled  by  the  geomagnetic  field. 
Although  the  movement  of  ionization  is  geomagnetically  controlled  at  all 
heights  above  150  km  (or  even  less),  it  is  convenient  to  set  the  lower 
boundary  of  the  magnetosphere  at  about  1500  km  where  the  neutral  particle 
density  becomes  predominant  over  the  ion  density.  The  magnetic  field, 
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Figure  10.  A  typical  daytime  midlatitude  electron-density  profile. 
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and  thus  the  magnetosphere,  end  at  .bout  ten  earth  radu  on  th  * 

°f  *"th  -  “  -*  o«  the  night  side  <see  ’’  ^ 
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appendix  c 


DERIVATION  OF  EQUATION  (I) 

The  development  given  here  roughly  parallels  that  given  In  Sections 

.3.1  and  5.3.2  of  Clemmo.  and  Dougherty  (1969).  Here,  however,  the 

effects  of  positive  ions  are  Included  and  some  different  notation  Is 
used . 

»e  consider  a  plasma  that  Is  unbounded,  cold,  and  colllslonless 

end  that  In  the  unperturbed  state  Is  uniform,  neutral  and  at  rest.  All 

wave-field  quantities  are  assumed  to  vary  as  exp  (iwt  -  ik-7)  where  w  ls 

the  angular  wave  frequency,  k  is  the  vector  wave  number  and  7  the  posi-' 

"ion  VtC£°r-  F°r  the  6leCtrons  -  «n  write  down  the  linearized  pertur- 
bational  equation  of  motion  as 

i,jmeVe  =  -e(E  +  v  x  5  } 

“  o 

Here  m^  ls  the  mass,  -e  is  the  charge,  v^  ls  the  perturbation  velocity 
E  is  the  average  electric  field  and  5^  the  static  magnetic  field  An  ’ 
analogous  equation  describes  the  ion  motion.  Ho.  letting  B  -  (0  0  B  ) 
and  noting  that  the  perturbation  current  is  J  c  He(7  v  ,  f°r  assumed! 
egual  unperturbed  number  densities  H  o,  electrons  and  ions,  ,h.  a  com- 
ponent  of  the  current  can  be  written  as 
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y  Jz 


1/2 


Here  Upe  =  (He2/eome)  ls 


the  electron  plasma  frequency  and  e>ce  =  e  Bo/me  is  the  electron  gyro 
requency.  Again  there  are  similar  definitions  for  „  ,nd  „  wlth 
only  the  mass  subscripts  changed.  P  Ci 
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where  4  and  eQ  are  the  free-space  permeability  and  permittivity  respec¬ 
tively,  and  e  is  defined  to  be  the  equivalent  relative-permittivity 
tensor  of  the  medium.  Substituting  for  j  and  solving  for  e  we  obtain 
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The  Maxwell  curl  equations  are  now 
V  x  E  =  -itoB 

and  V  x  B  =  i • 

Eliminating  5  and  performing  some  algebra  we  obtain  the  wave  equation 
k^E  -  k (k  •  E)  -  k^  £‘E  =  0 

where  kQ  =  to/c  is  the  free  space  wave  number.  Restricting  lc  to  lie  in 
the  x-z  plane  and  writing  out  this  last  equation  by  components  we  obtain 
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Dividing  through  by  -k  and  noting  that  k  ^  k  n  sin  0  and  k  = 

°  X  O  z 

ko  n  cos  0  ,  where  n  is  the  refractive  index  and  0  is  the  angle  between 
k  and  Bq,  this  equation  becomes 
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which  is  equation  (1)  of  the  text. 

Note  that  for  the  low-frequency  case  where  w  «  the  following 

approximations  can  be  made; 
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Here  c  is  the  vacuum  speed  of  light  and  va  is  the  Alfve'n  speed  defined  by 
2  2 

VA  =  ^o^o^mi  ’  **  ^as  136611  assumed  that  v^«  as  is  the  case  in  the  F 

region  of  the  ionosphere. 
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